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SONIC ENVIRONMENT OF AIRCRAFT STRUCTURE 
IMMERSED IN A SUPERSONIC JET FLOW STREAM 


Wiley A. Guinn, Frank J. Balena, and Jaak Soovere 
Lockheed-California Company 
Burbank, California 


SUMMARY 

Results of a study that was performed to make an assessment of the 
technological basis for using a small model to determine the sonic environ- 
ment on aircraft structure immersed in a supersonic jet flow stream is 
reported herein. 

Background information is given that pertains to noise source consider- 
ations, selection of test conditions, and resolution of hydrodynamic and 
acoustic pressure fields. 

Test, data requirements that are needed to make sonic fatigue, crack 
growth, interior noise and equipment analyses are defined. Data reporting 
formats are illustrated and data, applications are discussed. 

Instrumentation requirements for data acquisition, storage and reduction 
are given. Detailed discussions are given that pertain to pressure transducer 
and tape recorder requirements. Methods are given for computing cross- 
correlation coefficients from autocorrelation and space-time cross-correlation 
function plots. Instrumentation used for a SCAT 15-F sonic environment test 
is discussed and methods used for data reduction are explained. A list of 
methods for improving data acquisition and reduction for future sonic 
environment test programs is given. 

Methods are given for scaling of model data to full-size aircraft 
conditions. A literature search failed to provide a set of corresponding 
model and full-size engine test data associated with supersonic jet flows 
that was suitable for evaluating the scaling procedure. Therefore, subsonic 
jet test data for S-3A, L-1011, and a V/STOL configuration are compared. 

The SCAT 15-F model test data are scaled to dimensions and operating charac- 
teristics of a current SST duct burning turbofan engine concept. 

Appendixes include a derivation of the response spectral density equa- 
tion, a synopsis of the literature search, a listing of the computer program 
used for data scaling and computer program outputs for S-3A, L-1011, V/STOL, 
and SCAT 15-F test data. 



1 . INTRODUCTION 


A supersonic cruise aircraft research (SCAR) program was initiated by 
NASA in 1972 to develop technology for an advanced supersonic transport. 
Prediction of the sonic environment on aircraft surfaces that is caused by 
high velocity jet flows is one technology area that has been identified where 
advances are required. The accuracy of analytical prediction methods or 
scaling of model test data to full size aircraft dimensions has not been 
established. However, the sonic environment must be known in order to per- 
form analysis of sonic fatigue, crack growth, equipment environment and 
interior noise. 

One of the SCAR program aircraft concepts that requires sonic environ- 
ment evaluation is an over-the-wing engine configuration that gains lift 
for slow-speed flight by using Coanda turning of the jet stream. To achieve 
Coanda turning, the jet flow stream must be attached to the upper wing sur- 
face. The resulting thermal-acoustic environment on the parts of the wing 
surfaces that are immersed in the flow field will reduce the structural life 
of a typical wing structure. Therefore, special designs are required to with- 
stand the adverse environment. Design methods for skin-stringer-type struc- 
tures that will withstand the thermal-acoustic environment are given in 
References 1 through 3. However, more efficient designs may use hat-stiffened 
skins (Refs, k and 5), thermal tiles (Ref. 6) or ceramic composites. 

The current study was performed to investigate problems encountered in 
conducting model tests for supersonic jets and to evaluate accuracy of the 
test results. Background information leading to the study are given in the 
following paragraphs . 


1 . 1 Noise Source Considerations 

The degree to which model test data is comparable to full-size air- 
craft dimensions is dependent on several factors. These include the following 
noise source considerations: 

• Shock cell noise 

• Crackle 

• Internally generated noise 

•' Aircraft speed 

Shock cell noise is generated by supersonic jets when the exhaust flow 
is not fully expanded (Refs. 7 through 12). In cold jets a pure tone noise 
called screech is generated, whereas in hot jets the shock noise is broad 
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band (Ref. 13). Because, of difficulties associated with designing and 
manufacturing a practical convergent-divergent nozzle which operates at fully- 
expanded conditions during the entire flight mission, shock noise may occur 
during some portion of the flight (Refs. 9 through 12). Shock cell noise 
intensity and frequency are fvinctions of the nozzle pressure- ratio. ThereT- 
fore, as the pressure ratio changes with altitude, shock cell noise may sweep 
through the frequency range of predominant response. Consequently, premature 
structural failures may occur as a result of shock cell noise unless it is 
accoiinted for in the design phase (Ref. 9)* Model tests to investigate true 
effects of shock cell noise must be conducted with hot jets. 

Crackle may be significant in high-velocity jets (Ref. 13). When this 
phenomenon occurs, the noise signature has a distinctive bias toward, high- 
amplitude, positive, short-duration peaks. This results in increased ske-wnes 
of the noise probability density distribution with increasing jet velocity. 
Therefore, when crackle occurs , a substantial number of positive peaks may 
exceed the 3a rms peak level. 

Internally generated noise associated with rotating machinery, combus- 
tion, or high-velocity flows over obstructions is not simulated in a model 
test. If these noise sources are ■ significant , model test data may be 
misleading. 

The effect of aircraft speed on the sonic environment of a structure 
immersed in a jet flow stream is dependent on the relative importance of '■ 
hydrodynamic and acoustic-pressure fluctuations on the structure. If the 
acoustic pressure field is predominant, the structure environment is likely 
to decrease with aircraft velocity because noise generated by a jet is a func 
tion of the relative jet velocity (Vj. = Vj - Vg,). If the hydrodynamic pres- 
sure field is predominant, the aircraft velocity may not significantly alter 
the sonic environment since the jet velocity relative to the wing surface is 
virtually unchanged. Wind tunnel, sled, or whirling model tests are required 
to evaluate the effects of aircraft velocity on sonic environment 
characteristics. 


1.2 Selection of Test Conditions 

Before a meaningful, test can be conducted, aircraft operating conditions 
that are likely to establish noise design criteria must be determined. 

Static takeoff thrust generally produces the highest noise on a struct-ure 
in and adjacent to a jet flow stream. Thus, structure that is designed to 
withstand takeoff noise can usually withstand the sonic environment , for other 
operating modes. Reverse thrust also produces high noise levels that- may be 
predominant on some areas of the structure. Although jet noise decreases as 
aircraft velocity increases, the lower sonic environment may produce signifi- 
cant structiaral damage because of the relatively long length of exposure time 



during . flight . The relative importance of each test condition must he 
determined to ensure that test resiilts obtained will provide the correct 
environment for design of the structure. 


1. 3 Resolution of Hydrodynamic and Acoustic 
Pressure Fields 

Fluctuating pressures in a supersonic jet are composed of hydrodynamic 
and acoustic pressures. The impinging and attached jet flow surface pressure 
fluctuations (Refs. l4 and 15) and separated flow pressures (Refs. l6 through 
l8) which may occur over the trailing-edge control surfaces due to adverse 
pressure gradients differ from the acoustic field pressure in convection 
velocity and correlation signature. Therefore, it may be possible to resolve 
the two pressure fields by narrowband space-time correlation coefficient 
analyses if the local flow convection velocity differs from the speed of sound 
in the jet flow field. Generally, flow pressure fluctuations associated 
with boundary layers (Refs. 19 through 22) and with separated flow pressure 
exhibit a reduced coupling with the structure relative to that produced by the 
acoustic field. Therefore, unless the degree of coupling is considered when 
making structural design analyses, overdesign of the structure and correspond- 
ing weight increases are likely to occur. 


k 



2. TASK I - TEST DATA REQUIREMENTS 


Predictions for sonic fatigue, crack growth, equipment vibration and 
interior noise analyses require that response of the structure be determined. 
Frequency ranges normally investigated are: 

• Sonic fatigue and crack growth (50 to 1000 Hz) 

• Equipment vibration (50 to 2000 Hz) 

• Interior noise (50 to 10,000 Hz) 

The following sections define test data parameters, show data presenta- 
tion formats, and discuss data applications that are required to determine 
structural response caused by a random pressure loading on an aircraft 
structure . 


2 . 1 Data Requirements 

Noise Contours . -An OASPL distribution over the surface of an aircraft 
structure provides an indication of potential noise areas. The following 
guideline (based on rule-of-th\jmb estimates) can be used for making an assess 
ment of potential problems . 

• OASPL > 120 dB: Interior noise and equipment vibration problems may 

exist 

• OASPL > 150 dB: Sonic fatigue and crack-growth problems may exist 


Noise Spectra. -Spectral noise contours for the octave-band center fre- 
quencies over the range of 63 to 1000 Hz are generally sufficient to allow a 
designer to evaluate integrity and estimate weight of an aircraft structure. 
Octave-band, one-third-octave-band, or narrow-band frequencies provide suffi- 
cient information that can be used by empirical predictions to predict 
interior-noise and equipment-vibration environments. 


Correlation Functions and Spectral Densities . -Cross-spectral density of 
the excitation is required to compute response of a structure by the normal 
mode method. The cross-spectral density can be determined directly from a 
noise signal. However, when this is done, signal -phase relations are lost. 
Therefore, autocorrelation and cross-correlation functions are normally deter 
mined and Fourier Transforms are used to compute the spectral densities. 



Distribution Functions . -The probability density function indicates 
percentage of time that a random signal dwells between two amplitude limits. 
Two probability density distributions are commonly used in performing sonic 
fatigue analyses. They are: 

• Gaussian distribution . 

• Rayleigh distribution 

Jet noise is generally considered to have a Gaussian distribution that 
is defined by Equation (l). 


P(x) 


1 -x^/2o^ 

" ^ 


( 1 ) 


p(x) = probability density 

X 

2 
a 

a 

Instantaneous stress peaks of an aircraft structure that is subjected 
to a narrow-band random signal are generally considered to have a Rayleigh 
distribution that is defined by Equation (2). 


= instantaneous value of the noise signal with zero mean 
= variance =<x^(t)> 

/ 2 

= standard deviation = v<x (t)> 


p(x) 

X 


2 2 
X -X /2a 
= — e 

a 

=. instantaneous value of the envelope of noise peaks 


( 2 ) 


Assumptions are often made that excitation signal and stress response 
have Gaussian and Rayleigh distributions, respectively. Increased confidence 
in accuracy of fatigue analyses may be established by analyzing a random sig- 
nal to show that its distribution compares with the Gaussian and Rayleigh 
distributions . 


Model Data . -Atmospheric conditions, model geometry, flow conditions and 
pressure transducer characteristics should be reported. Also, flow-field 
boundaries, convection velocities and boundary-layer thickness may need to be 
determined. The use of these data will become evident in the subsequent 
discussions. 
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2. 2 Data Reporting Formats 


Table 1 gives a summary of the test data requirements that were defined 
in Section 2.1 and shows why each type of data is needed. The following 
paragraphs depict the format for data presentation. 


Noise Contours . -Figure 1 is a typical example of takeoff OASPL contours 
for an arrow wing supersonic transport. The values shown are for engines 
that are equipped with high-attenuation mechanical suppressors. The spectral 
level contours for the 63- to 1000-Hz frequency range octave-band center 
frequencies can be presented in a similar manner. 


Noise Spectra. -Figure 2 shows typical octave-band, one-third-octave- 
band, narrowband (20-Hz) and spectrum (l Hz) plots. The abscissa of these 
plots can be changed to Strouhal nimibers by application of Equation (3) for 
noise of a jet flow stream 


SN=^ (3) 

j 

SN = Strouhal number 

f = band center frequency 

D = diameter of nozzle 

V. = Jet velocity 
J 

The noise levels of one type of spectrum (e.g., one-third-octave band) 
can be converted to another type spectrum (e.g., octave bands) by application 
of Equation (4). 


SPL = SPL^ + 10 (4) 

r 

SPL = sound pressure level for desired bandwidth spectrum 

SPL^ = sound pressvure level for reference spectrum 

Af = desired frequency bandwidth 

Af = reference bandwidth 

r 
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TABLE 1. UTILIZATION OF SONIC ENVIRONMENT TEST DATA 
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Figure 1 . Typical Engine Noise Contours 
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Figure 2. Typical Noise Spectra Plots 
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Correlation Density Functions and Spectral Densities . -Figure 3 shows a 
typical power spectral plot, this plot can be converted to sound pressure 
levels for any bandwidth by Equation (5)- 


SPL = 10 log 


10 — * 10 logio Lt 


p = ordinate value of Figure 3 


( 5 ) 


P 


2 

r 


= reference pressure squared 


Figure h shows typical autocorrelation and cross-correlation plots. 
Section 3.2 contains equations for computing correlation coefficients from 
these curves. Section 2.3 describes methods for determining the spectral 
density by using the correlation coefficient values. 


Distribution Functions .-Figure 5 gives plots of the probability density 
Gaussian and Rayleigh distributions that are defined by Equations (l) and (2), 
respectively. The curves shown are normalized with respect to the standard 
deviation. 


2.3 Utilization of Data 

Sonic Fatigue . -The first step in sonic fatigue analysis is to establish 
the design life of a structure at the highest noise level. This is achieved 
by studying the aircraft utilization pattern. Parameters included are take- 
off, landing, flight profiles, ground taxi, and static aircraft test noise 
levels. A general procedure used for sonic fatigue analyses (Ref. 9) is to 
account for noise-reduction levels during the ground run and flight by com- 
puting the equivalent damage duration at static takeoff noise levels. This 
procedure assumes that the nature of the flight and takeoff acoustic environ- 
ments remain essentially the same. 

It is current practice to require sonic fatigue proof testing of any 
novel structures, such as those required for high thermal-acoustic environ- 
ments, because structural details cannot be accounted for by analyses. These 
tests are usually conducted in an acoustic progressive wave test facility 
that provides for adjustment of spectrum levels and shapes for a predefined 
correlation function. 

Two basic methods are used for performing sonic fatigue analyses. The 
first is an empirical method that is based on design charts. The other 
method is the normal mode approach based on the procedure given in Ref- 
erence 23. Current sonic fatigue analyses are based on a predominant 
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single-mode response. The analyses are used in conjunction with random 
fatigue data for critical joints. Most random fatigue data are based on the 
Rayleigh distribution of stress peaks which is obtained from a single modal 
response to a broadband Gaussian- type random excitation. A reasonable 
approximation to the Rayleigh distribution occurs when the multimodal response 
falls within a frequency band that has an upper limit approximately twice 
that of the lower limit (Ref. 2h) . The current trends in generation of ran- 
dom S/N data by coupon testing is to use a broader spectrum of excitation 
to include the contribution from the higher modes. 

The empirical analysis method of obtaining data is expensive and is 
usually restricted to simple rectangular structiures such as skin-stiffener 
structure, simple honeycomb panels and beaded panels. After the panel dimen- 
sions are selected, the panel frequency is computed. The structural life is 
determined by assuming the single-mode response and using the equivalent 
damage duration for static takeoff levels. This process is repeated until a 
suitable design has been achieved. 

The basic normal mode approach equation for the response spectral density 
at a point (x,y) is given by Equation (6). 


G(x,y;oj) 


EE 

v»= n o ?= T 


1 1. 


2 - 2 2 2 
M (o) -oj +2i6 0 ) to) M (u) -co +2i6 to to) 
rr rr ss ss 


( 6 ). 


X f^(x,y)fg(x,y) 



•^r ^ ^1 ’^1 ^ ^s ^ ^2 ’^2 ^ ^ ^ 


G(x,y;to) 


fj.(x,y) ,fg(x,y) 

f(x^,y^) 

f(x2>y2) 

Gp(S,h;o)) 


M ,M 
r ’ s 


to 

r 


jto 


s 


response spectral density 

normal mode deflection at point (x,y) for modes r and s 
normal mode deflection at point 
normal deflection at point 

excitation cross-spectral density at separation dis- 
tances 5 and n between points x^jy^) 

generalized mass corresponding to r and s modes 

natural frequency corresponding to mode shapes r and 
s, respectively 
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“ = forced frequency 

6^)6g = viscous damping factor corresponding to modes r and s 

Derivation of Equation (6) is given in Appendix A. The first two factors 
following the siunmation signs are the receptances of the system for modes r 
and s, respectively. The double integral term represents the degree of 
coupling between the excitation and structural response. The equation is 
applicable to jet noise, turbulent boundary layer excitation (Ref. 20), and 
separated flow excitation (Ref. 17)- For turbulent boundary-layer excitation, 
the cross-spectral density can be computed by Equation (?)• 


-iw5/U 

Gp(C,n;o3) = Gp(cj) |pp(C,0,T;m) ( |pp(0,n,T;o)) I e 


( 7 ) 


direct spectral density 

longitudinal narrowband space-time correlation coefficient 

lateral narrowband space-time correlation coefficient 

subscript designating excitation quantities 

longitudinal separation distance between transducers 
(e.g. transducers 1 and 2) 

n = lateral separation distance between transducers 

= convection velocity of the flow stream 

Narrow-band correlation coefficients for a traveling acoustic wave at grazing 
incidence can be computed by use of Equation (8) or from test data (Sec- 
tion 3.2). 


Gp(m) 

Pp(C,0,T;w) = 
Pp(0,n,T;u) = 

P 

C 


p(5, n, T, w) = cos (t - ?/a) 


( 8 ) 


a = speed of sound in the flow field 



Considerable simplification of Equation- (6) is obtained by making 
assumptions that are commonly used when making sonic fatigue analyses. The 
assumptions are: 

• A predominant single-mode response 

• A fully correlated excitation across the panel 

• A constant excitation spectrum level. 

The double area integral in Equation (6) is reduced to an integral of mode 
shapes if the pressure field is assumed to be fully correlated over the 
panel area. This assumption is correct for normal incident acoustic waves 
and results in very small errors for a fundamental mode progressive wave. 
Therefore, the mean-square response obtained by integrating Equation (6) 
with respect to circular frequency (m) results in Equation (9a). 


<y (x,y,t )> = 


_ 

i+M 03 6 

r r 


(9a) 


A similar expression to Equation (9a) was developed by Miles (Ref. 25). The 
expression (Equation 9^) is based on the static displacement yo(x,y) at 
point (x,y) due to a unit pressure on the structure. 


2 TT 2 

<y ( x,y,t)>= m^Gp(o))y^(x,y) 


(9b) 


The corresponding’ expression for mean-square panel stress is given by Equa- 
tion (lO). 


2 TT 2 

<a ( x,y,t)> = u)^Gp(o))o^(x,y) 


( 10 ) 


a = static stress at point (x,y) on the structure due to a \init pres- 
sure over the structure 

The panel stress is used with random fatigue data for representative struc- 
tiire to determine fatigue life of the structure (Ref. 26). 
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Crack Growth . -Crack-growth analyses (Ref. 27) are "based on a modified 
Rayleigh Ritz method of assumed cracked-panel modes. Initially the panel 
response spectral density is computed "by Equation (6). Expressions for com- 
puting stress spectra are then developed by using the assumed cracked-panel 
modes. Baseline panel crack-growth data due to a random loading are obtained 
from electromagnetic shaker-excited coupon specimens. The baseline crack- 
growth test data are used in conjunction with computed stress spectra to 
predict crack growth. 

Equipment Vibration . -Dynamic characteristics of structural vibration 
reflects the combined effects of sonic environment and structural response 
characteristics. Complexity of the aircraft structure makes a theoretical 
prediction method impractical for making engineering analyses. Therefore, 
empirical methods are used. These methods are based on the use of correla- 
tion curves. The curves are established by correlating acceleration response 
levels that are measured on primary structure of existing aircraft with the 
aircraft sonic environment (Refs. 28 and 29)- The problem approach is: 

• Division of the aircraft into zones of approximately equal-vibration 
response levels 

• Estimation of the octave-band acoustic levels over the aircraft flight 
conditions of interest 

• Use of response correlation cxurves (Refs. 28 and 29) 

• Prediction of frequency-dependent vibration spectra 

Typical vibration zones for a jet-powered subsonic airplane are shown 
in Figure 6. Figure 7 illustrates the procedure for converting the sonic 
environment of each zone to acceleration spectral density. A typical environ- 
ment for the outboard wing area is shown in Figure 8 (Ref. 29). The spectral 
density levels are used as standards for equipment qualification test levels. 


Interior Noise . -Interior noise in passenger-occupied areas of an aircraft 
that is associated with jet noise is maximum at takeoff. It is still present 
during flight at a level comparable to turbulent boundary-layer noise for the 
lower-frequency region (Ref. 30). Once the exterior noise levels on the 
fuselage have been determined, interior noise analysis includes obtaining 
values for the following quantities. 

• Transmission loss of the structure and acoustic treatments 

• Interior absorptivity 

• Interior equipment noise 
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Figure 6. Vibration Zones and Their Location 
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• Mechanical vibration noise generation 

• Size and shape of interior 

Computation of the interior noise is an intricate task. Currently, 
analytical prediction methods have not been developed that can be used to 
accurately predict the interior noise. Methods that are representative of 
the technology are given in References 22 and 31- Empirical methods for 
interior noise predictions rely heavily on test data banks that have been 
compiled which give noise levels in existing aircraft along with the corre- 
sponding acoustic treatments. The aircraft data are supplemented by mounting 
a representative section of the aircraft fuselage (e.g., structure plus 
acoustic treatment plus interior trim) between two reverberation rooms and 
measuring the transmission loss. Normally, the noise on the exterior side of 
the panel has approximately the same spectral content as the jet or turbulent 
boundary- layer noise, but it is a normal incident wave which does not simulate 
the degree of coupling for a turbulent flow field. Nevertheless, comparison 
of different acoustically treated panel configurations provides a relative 
comparison of the panel noise reduction characteristics. 
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3. TASK II - DATA ACQUISITION AND REDUCTION 


3 . 1 Instrumentation and Procedures 

Pressure transducers and tape recorders that are used to record the noise 
for a Jet model test must be carefully selected. Required characteristics for 
each of these are discussed below. 


Pressure Transducers .-Selection of pressure transducers for a sonic 
environment test must include the following considerations. 

• Environmental conditions to which the transducers are exposed (e.g. 
temperature, humidity, and vibration) 

• Size of the sensing element, 

•. Dynamic range 

0 Frequency response. 

Temperature Environment: The temperature environment in a hot Jet flow 
stream is a formidable requirement for pressure transducers. Pressure trans- 
ducer manufacturers have developed several transducers for measuring pressure 
fluctuations in a high temperature Jet flow stream. The suitability of these 
for use on hot Jet model tests .remains to be determined. Limitations for 
various types of the high temperature pressure transducers include the 
following : 

• They cannot be flush mounted 

• They require water cooling 

• They have insufficient frequency response 

The acc\aracy of these transducers needs to be determined. The rationale 
for this statement is based on a comparison of measurements made by 12 differ- 
ent low-temperature pressure transducers (Ref. 32) in a wind tunnel. Data 
recorded by the various transducers for Mach numbers of 1.6 to 2.5 showed 
significant differences. A similar test has not been conducted for high- 
temperature transducers, but it is anticipated that a comparison of data 
recorded by different models would result in large discrepancies. 


Size of Sensing Element: Finite size of a transducer sensing element 
limits its space resolution of a pressure field (Ref. 33). As the value of 
the quantity coR/U^ increases, there is a corresponding increase in measurement 
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error (Figure 9)- For a given jet flow stream, the circular frequency (03) 
and the convection velocity (U^) are fixed. Consequently, space resolution 
can he improved only by making the transducer sensing element radius (R) 
smaller. 


Dynamic Range: The dynamic range of a pressure transducer must be com- 
patible with the magnitude of pressure fluctuations that are to be measured. 
The lower level of the range is limited by the signal-to-noise ratio and the 
upper level is limited by clipping of the signal. 


Frequency Response: The transducer frequency response required for model 
testing depends on the model scale (Section it.l). As the model size, is 
decreased, the frequency range to be measured increases. Transducers that 
are suitable for measurement of high-frequency noise need small sensing ele- 
ments to ensure good frequency response. However, frequency response and 
sensitivity of a transducer vary inversely. Therefore, the most suitable 
transducer for making sonic environment measurements is the one with the 
smallest sensing . element that has sufficient sensitivity for the intensity 
levels being measured. 


Data Storage . -Test data are generally stored on magnetic tape. This can 
be accomplished by recording data in the direct or the FM mode. The mode to 
be used depends on the frequency bandwidth to be measured and the manner in 
which the data will be analyzed. 


Direct Recording Mode: The direct mode permits measurements up to 6 OO kHz 
in the intermediate band mode of operation and to 2 MHz in the wideband mode 
of operation. Disadvantages of the direct mode are poor low-frequency 
response, complexity of frequency response corrections for time expansion, 
amplitude instability (commonly referred to as dropout) at very high fre- 
quencies and low signal-to-noise ratio. The poor low-frequency response will 
not be a problem if the model scale is sufficiently small so that measurement 
of frequencies below approximately 1*00 Hz are not required. Time expansion 
is not required if a spectral analyzer is used for data processing that has a 
sufficiently wide bandwidth so that data can be reproduced at the same speed 
at which it is recorded. Amplitude instability can be minimized by using 
high-quality magnetic tape and keeping recorder heads , guides and other parts 
of the recorder that come in contact with the tape scrupulously clean. The 
low signal-to-noise ratio is a definite limitation. 


FM Recording Mode: The FM mode has good amplitude stability (virtually 
insensitive to dropouts) and low-frequency measurement capability (down 
to dc). Wideband Group 2 FM recording permits measurements of frequencies 
ranging from dc to 500 kHz. Time expansion can be accomplished by recording 
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Resolution of Frequency Spectral Density for a Round Transducer 




at a high tape speed and playing back at a low tape speed. When this is 
accomplished in the FM mode, minimal frequency response corrections are 
required in comparison to those for the direct mode of operation. This time- 
expansion capability is useful when analyzing transient signals or when the 
measured data bandwidth is wider than that of the data-reduction analyzer. 
Group I FM recordings have a singal-to-noise ratio that is approximately 15 dB 
greater than the direct mode and the Group II signal-to-noise ratio is approx- 
imately the same as that of the direct mode. The Group I mode can be used if 
frequencies of the data to be measured do not exceed approximately 80 kHz. 

The FM mode of operation is considered to he most favorable in light of 
the aforementioned considerations. 


Phase Calibration . -When cross-correlation plots are to be made, a phase 
calibration of all data channels that are to be used for cross correlations 
is necessary. If this calibration is not performed, the cross-correlation 
functions will include initial phase differences which result in a time shift 
of the entire function. It is good practice to record all data to be corre- 
lated on either even- or odd-numbered tape recorder channels . This eliminates 
the possibility of errors caused by differences in the recording head loca- 
tions for even and odd channels.. 


Cross-Correlation Coefficient . -Figure 10 illustrates the manner in which 
narrowband cross correlation coefficients are determined. The autocorrelation 
functions are determined by taking the time average of the product p(t)p(t+x) 
where x is the delay time. Rp^(x) and Rp^Cx) are typical narrowband plots for 
pressure measurements at locations Pp and Pp, respectively. The cross- 
correlation plot is determined by taking the time average of the product 
Pl(t)p2(t+x) . This plot is the lower plot in Figure 10. The cross-correlation 
coefficient is obtained by dividing the peak amplitude of the cross-correlation 
plot corresponding to the delay time xp p "the square root of the product 
of the amplitudes of the autocorrelation functions for x=0. Figure 10 shows 
the narrowband autocorrelation functions and the cross-correlation function to 
be slowly decaying periodic functions. However, correlation plots for wide- 
band random signals decay rapidly to zero as the value of x increases. 


3-2 SCAT 15-F Model Test Data Analyses 

A schematic diagram of the data acquisition and data-processing system 
used for the SCAT 15-F model test (Ref. ) are shown in Figure 11. 


Pressure Transducers .-Three different models of pressure transducers were 
used for measuring the sonic environment of the. upper wing surface during the 
SCAT 15-F model test. They were Bruel and Kjaer (B&K) 4138 microphones, 

Kulite VQL-250-25 transducers and Piezatronics 112A02 pressure transducers. 
Characteristics of these transducers are listed in Table 2. 
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Figure 11 . Data Reduction and Analysis System 
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(TABLE 2. SCAT 15-F iTEST PRESSURE ' TRANSDUCER CHARACTERISTICS 



B&K 

4138 

KULITE 

VQL-250-25 

PIEZATRONICS 

112A02 

Diameter 

3.1T5xlO"8m 

6 . 35xIO“8m ■ 

5.537xI0“8m 

Dynamic Range 

7 6-1 68 dB 

— 

131 to 211 

Frequency Response 

7-i 40 kHz 

— 

— 

Resonant Frequency 

— 

35 kHz 

250 kHz 

Vibration Sensitivity 

IG = 80 dB 

IG = 100 dB 

0.002 N/m2/G 

Thermal Sensitivity 

0.0028 dB/°C 

0.011% FS/°C 

0.011% FS/°C 

Static Pressure Sensit. 

-1 dB/ATM 

NA 

NA 


Initially three B&K microphones were flush mounted in the wing surface 
of the model (Figure 12).. However, as the jet velocity was increased during 
the first test condition when the nozzle was located at position 1, the 
diaphragm of one of the microphones was destroyed. Since time allotted for 
conducting the test was two weeks and the primary objective of the test pro- 
gram was to determine far-field noise reductions that can be attained 
through shielding of a jet noise source by an arrow wing structure, the B&K 
microphones had to be replaced by transducers that were readily available. 

Two Piezatronics transducers and one Kulite transducer appeared to be the 
best that were available. Therefore, they were flush mounted in the wing sur- 
face. These transducers are extremely rugged. The Kulite transducer is a 
solid state sensor that is rated for l.T24xlo5 N/m^ with a maximum usable 
pressure of 3.44 Tx 105 N/m^ and the Piezatronics pressure transducers can 
withstand a maximum static pressure of 1.3x10^ N/m^. 


Environmental Conditions . -Environmental conditions did not appear to 
have a significant influence on the choice of transducers. Since the 
test was conducted in an anechoic room, humidity was not considered to be a 
problem. Model weight and rigidity of the model support were believed to be 
sufficient to prevent excessive vibration levels that would affect the noise 
measurements. Air supply to the nozzle was near ambient conditions. There- 
fore, the jet temperature was not considered to affect the transducer sensi- 
tivities. However, the fully expanded jet static temperature was about 
-157 degrees Celsius and may have been a factor that contributed to 
failure of the B&K microphone. 
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PRESSURE TRANSDUCERS 


PIEZATRONICS 112A02 | 
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Figure 12. SCAT 15-F Sonic Environment Test Configuration 
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Finite Size of Sensing Element . -The sensing element for the Kulite and 
for the Piezatronics transducers was 6.35x10-3 and 5.53TxlO-3m, respectively. 
Therefore, if the convection velocity (U^) is considered to he 0.62 of the jet 
velocity, convection velocities for the Mach 2.5 (Vj = 550 m/s) and 
Mach 1.5 (Vj = 427 m/s) are 340 and 265 ni/s, respectively. Correction 
values (10 log^Q for the 80 kHz upper frequency are: 


u)R/U^ 

(Fig- 9) 


Mach 1.5 Nozzle 
Kulite Piezatronics 

6.05 ' 5-27 

0.0067 0.0103 


Mach 2.5 Nozzle 
Kulite Piezatronics 

4.70 4.1' 

0.0124 0.0136 


10 log. 


<t>m 


10 ' 


-21. Jh dB 


- 19.87 dB 


- 19.06 dB 


- 19.87 dB 


As can he seen, the finite size effect for the transducers used for the 
SCAT 15-F test is large. 


Dynamic Range. -The sonic environment on the wing surface of the SCAT 15-F 
model was estimated to he within a range of 100 to 160 dB. Table 2 shows that 
the dynamic range of the B&K 4138 microphone is 76 to I 68 dB. However, the 
sonic environment may have exceeded the upper limit of the dynamic range and 
contributed to failure of the microphone. The dynamic range of the Pieza- 
tronics and Kulite transducers is suitable' for the higher intensity environ- 
ment encountered' in the test. 

Sensitivity of the Piezatronics transducers used for the SCAT 15-F model 
test were considered to be marginal for the range of pressures measured. 

Table 2 shows the lower limit of the dynamic range to be 131 dB. Therefore, ' 
internally generated noise of the measuring system may have affected the lower 
intensity noise level measurements. Frequency response calibrations were not 
available for either the Piezatronics or Kulite transducers and means were not 
available for ' performing them. Therefore, the response was assiunedto be 
uniform with frequency. . 


Frequency Response. -Model scale for the SCAT 15-F model test was con- 
sidered to be 0 . 03 . Therefore, if the model jet velocity is considered to be 
equal to full-size engine jet velocity, fa, = 0.03 x see Section 4.1) . 
Measured noise levels covered the frequency range- from 50 to 80,000 Hz. Con- 
sequently, the corresponding full-scale frequency range was from I .5 to 
2400 Hz. It should be noted the the frequency range of interest for struct- 
ural analyses (50 to 2000 Hz) is well within limits of the measured noise 
levels . 
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Data Storage . -The FM mode was used for recording the SCAT 15-F model 
test data. Data were recorded at 3.048 m/sec tape speed on a 432 kHz carrier. 
Time expansion was accomplished hy playing back the tape at 0.38 m/sec on a 
54 kHz carrier. Therefore, the 80-kHz frequency was reduced to 10 kHz, and 
it was possible to reduce the data with a spectral analyzer that had a 10-kHz 
upper frequency limit. 


Data Reduction . -Upper-wing surface pressure data were recorded during the 
SCAT 15-F test runs 46P and 4TP at Mach 2.5 and runs POP, 31P, 32P, 33P, 36P, 
3TP, and 92P at Mach 1.5- The test parameters for each of the conditions are 
given in Table 3. The locations of microphones 10, 11, and 12 relative to 
the several jet locations used are shown in Figure 12. The initial data- 
reduction included one-third-octave-band analyses from 50 Hz to 80 kHz (Fig- 
ure 13) and narrowband analyses from 50 to 80 kHz (Figure l4). Later the 
narrowband data were plotted for 50 to l6 kHz (Figure 15) in order to better 
resolve the frequency content. 

In the process of reducing the data, a difference of 7-2 dB was noted 
between the pre- and post-calibration of the Kulite transducer (Location lO). 
These calibrations were recorded several days apart, and it was not readily 
apparent when the shift occurred. Calibrations were performed each day and 
used to verify the operation of each microphone system prior to each day of 
testing. However, the calibrations were not recorded on magnetic tape each 


TABLE 3. SUMMARY OF TEST CONDITIONS 


Run 

No, 

Nozzle 

Location 

Pressure 

Ratio 

Nozzle 

Mach No. 
\ 

Exit Velocity 
Meters/Second 

31P 

1 

3.67 

.1-5 

427 

32P 

2 

3.67 

1.5 

427 

33P 

3 

3.67 

1.5 

427 

2 OP 

4 

3.67 

1-5 

427 

■ 92P 

5 

3.67 

1.5 

427 

36P 

6 

3-67 

1.5 

427 

3TP 

7 

3.67 

1.5 

427 

46p 

6 

1.70 

2.5 

550 

47P 

7 

1.70 

2.5 

550 
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Figure 15 . Typical Narrow Band Spectrum (50 to l6,000 Hz) 




lay, in the interest of completing all of the scheduled tests within the 
allotted time period. The change in calibration level was not noted \mtil 
the post-calibration was recorded at completion of the test program. The 
calibration problem was investigated for the Mach 1.5 tests by plotting 
sound pressure level contours (Figure l6). First, the precalibration value 
was used and then the post-calibration value, was used. The contours which 
assumed the precalibration value to be correct were disjointed, whereas, 
those which assumed the post-calibration value to be correct was smoother. 

The contours given in Figure l6 were constructed from the transducer data 
grid shown- in Figure 17. 

The narrow-band plots such as the one shown in Figure 15 were used to 
determine frequencies for making correlation plots. The frequencies chosen 
are given in Table h. Typical autocorrelation and cross-correlation plots 
are shown in Figures l8 through 20. The plots are for run 31P at a 1000-Hz 
frequency. Figures l8 and 19 are the autocorrelation plots for transducers 
11 and 12 respectively. Figure 20 is the cross-correlation plot for 
transducers 11 and 12. The ordinate of plots l8 through 20 are given in 
terms of linear dimensions that are proportional to pressure squared. Since 
the autocorrelation and cross-correlation curves are used only for computing 
normalized cross-correlation coefficients, the conversion factor for convert- 
ing linear dimensions to pressure squared cancel out. Computation of the 
cross-correlation coefficient at the top of Figure 20 was accomplished by 
using the method illustrated in Figure 10. The narrow-band cross-correlation 
coefficients were determined for all of the frequencies defined in Table h in 
a similar manner. Tables 5 and 6 give the correlation summaries for the 
Mach 1.5 and 2.5 nozzles', respectively. The maximum correlation coefficient 
values given in column 5 are the pp(0,r) ,t ,o)) values that are used in 
Equation (T) to compute excitation cross-spectral density. 


3. 3 Methods for Improving Data Acquisition and Reduction 

Improvements listed below are with reference to the SCAT 15-F model test 
program. 

• A hot jet will better simulate a SST-type engine. Current SST 
engine concepts have exhaust velocities on the order of 823 m/s. 

• The pressure transducers should have smaller sensing elements to 
improve high-frequency space resolutions. 

• The transducer array should include longitudinal and lateral pos- 
itions so that true convection velocities can be determined. 
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Figure 17 . Contour Data Grid 
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TABLE U. CORRELATION MATRIX 


Run No. 

Transducer 

A 

Transducer 

B 

Correlation Frequencies 

- kHz 

31P 

11 

12 

1.0 

2.1 

1.8 


32P 

11 

12 

1.0 

2.1 

1.8 


33P 

12 

10 

1.0 

2. 1 

1.8 


33P 

. 12 

11 

1.0 

2.1 

1.8 


20P 

11 

12 

1.0 

2.1 

1.8 


92P 

12 

11 

1.0 

2.1 

6 . 6 


92P 

12 

10 

1.0 

2.1 

6.6 


36P 

11 

12 

1.5 

■ 2.1 

6.6 


37P 

12 

11 

1.0 

2.1 

6.6 

15 ' 

37P 

12 

10 

1.0 

2.1 

6. 6 

15 

37P 

11 

10 

1.0 

2.1 



46p 

11 

12 

1.0 

6.6 

9.0 

15 

47P 

10' 

12 

1.6 

6.6 

7.2 

15 

hlP 

11 

12 

1.6 

6.6 

7.2 

15 ■ 
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Figure l8. Autocorrelation Plot for Run 31P Transducer 11 







31P 1kHz Ri2, WO. 34 
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Figure 19- Autocorrelation Plot - for Run 31P Transducer 12 





COMPUTATION OF NORMALIZED 
CROSS CORRELATION COEFFICIENT 



Figure 20. Cross-Correlation Plot for Run 31P Transducer ll and 12 
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TABLE 5. MACH 1.5 NOZZLE CORRELATION SUMMARY 


1 



Maximum 







Correl. 

Correlation 

Maximum 

Peak 

Run 


Freq 

Function 

Correlation 


T 


No. 

Transducers 

(kHz ) 

(m) 


Coefficient 

( 

sec ) 

31P 

11-11 

1 

6.56 

X 

10 2 


- ■ 

0 



12-12 


7.33 

X 

10'2 



0 

10 ^ 


11-12 


2.98 

X 

10“2 

0.43 

1.46 

X 


11-11 

2.U 

2.71 

X 

l°l2 



0 



12-12 


3.70 

X 

10 ^ ! 



0 



11-12 


7.49 

X 

10-3 

0.24 

1.775 x 10 ^ 


11-11 

12-12 

1+.8 

4.11 

4.57 

X 

X 

1°I2 
10 ; 



0 

0 



11-12 


6.60 

X 

10 

0.15 

1.96 

X 

10,^ 

32P 

11-11 

1 

7.32 

X 

10~2 



0 



12-12 


6.86 

X 




0 



11-12 


3.86 

X 

10-2 

0.55 

1.30 

X 

10"^ 


11-11 

2 .k 

3.33 

X 

10"2 



0 



12-12 


4.85 

X 

10~2 ■ 



0 

10"^ 


11-12 


1.19 

X 

10-2 

0.30 

'O.65 

X 


11-11 

4.8 

5.51 

X 

10~2 : 



0 



12-12 


8.28 

X 

10-2 



0 



11-12 


1.22 

X 

10-2 

0.18 

0.80 

X 

10-^ 






-2 





33P 

10-10 

1 

: 3.10 

X 

l°-2 



0 



11-11 


'7.37 

X 




0 



12-12 
10-12 
11-12 ■ 


3.38 

3.3 

'7.37 

X 

X 

X 

? 

10:3 
10-3 
10 ^ 

0.10 

0.15 

4.2 

1.9 

0 

X 

X 

1°I4 

10 


10-10 

2.4 

8.08 

X 

loi 



0 



11-11 


8.64 

X 




0 


• 

12-12 


7.54 

X 

10 2 



0 

), 


10-12 


6.86 

X 

io:i 

0.09 

0.65 

X 

1°I4 


11-12 


2.21 

X 

10 

0.27 

0.45 

X 

10 


10-10 

11-11 

4.8 

5.46 

6.99 

X 

X 

1°I2 
10 , 



0 

0 



12-12 


7.87 

X 

'-°-3 
10 ^ 



0 

1, 


10-12 


.7.62 

X 

0.12 

1.9 

X 

10 1 


11-12 


2.54 

X 

0.03 

1.4 

X 

10-^ 


h3 
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Maximum 

Correlation 

Coefficient 


X 

( sec ) 


11-11 

12-12 

11-12 

1 

6.85 X 10 ^ 

6.46 X 10 ^ : 

2.24 X 10 

11-11 

12-12 

11-12 

2.4 

7.94 X 10"^ 
7.68 X 10 ‘Z 
2.91 X 10 ^ 

11-11 

12-12 

11-12 

4.8 

2.92 X 10“^ 

3.49 X 10 ^ ! 

1.27 X 10 i 

10-10 

11-11 

12-12 

10-12 

11-12 

1 

3. '87 X 10"^ 

7.94 X 10 ^ ! 

6.20 X 10~^ 

1.07 X 10 ^ 

2.4l X 10~ 

10-10 

11-11 

12-12 

10-12 

11-12 

■ 2.4 

i 

9.02 X 10”^ 

2.54 X 10 ^ 

5.59 X 10 ^ : 

7.11 X 10 :r 
. 4.32 X 10 ^ 

10-10 

11-11 

12-12 

10-12 

11-12 

6 . 6 

3.86 X 10"^ 

5.64 X 10 ^ 

; 5.08 X 10 ^ 

2.79 X 10 

9.65 X 10 ^ 

1 

11-11 

12-12 

11-12 

1.5 

4.38 X 10"^' 
' 4.85 X 10 

, 1.47 X 10 

11-11 

12-12 

11-12 

2.4 

’ 5.84 X I0~'i 

j 5.91 X 10 

9.65 X 10 ^ 

11-11 

12-12 

11-12 

6.6 

'4.94 X lO"^- 
2.71 X 10 ‘i. 

1 8.89 X 10 


1.T8 X 10 
1.83 X 10 


6.3 X 10' 


2.82 X 10 


4.77 X 10' 


44 
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TABLE 5. MACH 1.5 NOZZLE CORRELATION SUMMARY - Concluded 





Maximum 







Correl. 

Correlation 

Maximum 

Peak 

Run 


Freq 

Function 

Correlation 


T 


No. 

Transducers 

(kHz) 

(m) 


Coefficient 

( 

sec ) 


10-10 




-2 





37P 

1.5 

7.25 

X 

10 ^ 



0 



11-11 


. 6.54 

X 

10 p 



0 



12-12 


i 7.87 

X 




0 



10-12 


' 5.59 

X 

10 1 

0.074 


0 

0 


11-12 


4.42 

X 

10“^ 

0.616 


0 

0 


10-10 

2.4 

’3.77 

X 

10~2 



0 



11-11 

12-12 


: '5.65 

1 ^.95 

X 

X 

10 ^ 



0 

0 



10-12 


.2.10 

X 

10 2 

0.485 

1.66 

X 

10 1 


11-12 


2.92 

X 

10-2 

0.719 

1.65 

X 

10-^ 


10-10 

6.6 

7.68 

X 

10~2 



0 



11-11 


6.44 

X 

10-2 



0 



12-12 


2.07 

X 

10~2 






10-12 


1 2.79 

X 

10'^ 

0.07 

6.59 

X 

10 ^ 


11-12 


3.81 

X 

10"^ 

0.104 

8.. 63 

X 


10-10 

15 

i :5.84 

X 

10"2 



0 



11-11 


! 7.24 

X 


- 


0 



12-12 


■ 4.76 

X 

10-2 



0 



10-12 


1.27 

X 

10-3 

0.02 


- 



11-12 


2.5^ 

X 

10-3 

0.04 


- 


Repeat 



8.89 


P 





3TP 

10-10 

1.5 

X 




0 


11-11 


.4.83 

X 

10-2 



0 



10-11 


3.30 

X 

10-2 

0.54 

6.12 

X 

10 ^ 


10-10 

11-11 

2.4 

3.45 

4.72 

X 

X 

^°I2 

10 2 



0 

0 

), 


10-11 


1.52 

X 

10 2 

0.38 

1.3 

X 

10. 






•> 
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TABLE 6. MACH 2.5 NOZZLE CORRELATION SUMMARY 





Maximiun 






Correl. 

Correlation 

Maximum 

Peak 

Run 


Freq 

Function 

Correlation 


T 

No. 

Transducers 

(kHz) 

(m) 


Coefficient 

( 

sec ) 

46p 

11-11 

1 

8.00 

X 

10 J 



0 


12-12 


4.60 

X 





11-12 


4.65 

X 

10-2 

0.77 

2.5 

X 10 ^ 


11-11 

6 . 6 

7.11 

X 

l°l2 





12-12 


7.70 

X 

10 p' ; 



0 

1 


11-12 


: '4.83 

X 

10-2 

0.65 

.9 


11-11 

9 

4.34 

X 






12-12 

11-12 


3.05 

1.27 

X 

X 

10 ; 

10 ; 

0.35 

1.02 

-4 

X 10 






^-2' 





11-11 

15 

T.TT 

X 

10 p 





12-12 


4.57 

X 

10 



X 10-^ 


11-12 


6.25 

X 

0.11 

0.6 

4tp 

10-10 

1.6 

4.11 

X 

10-2 





11-11 


1.85 

X 






12-12 


8.81 

X 

lO'p! 



), 


10-12 


2.41 

X 


o.4o 

2.3 

X 1Q-_1 


11-12 


2.11 

X 

10 

0.52 

1.2 

X 10 


10-10 

6. 6 

' 1.80 

X 

10-2 
10 - 





11-11 


4.11 

X 





12-12 


3.84 

X 




-4 

X 10 1 


10-12 


6.35 

X 

10 ^ 

0.24 

io.6 


11-12 

• 


3.30 

X 

10 ^ 

0.08 

■0.5 

X 10 ^ 






-2 





10-10 

7.2 

■2.77 

X 






11-11 


6 . 60 

X 

10 2 





12-12 

10-12 


: 3.48 

3.3b 

X 

X 

10 i 

0.11 

3.3 

-4 

X 10 


11-12 


2.46 

X 

10-2 

0.51 

5.6 






-2 





10-10 

15 

3.33 

X 

10 





11-11 


2.36 

X 

10I 





12-12 


3.89 

X 

^ 0:3 

10 1 



-4 
X 10 , 


10-12 


5.08 

X 

0.l4 

0.4 


11-12 


1.09 

X 

10'^ 

0.36 

0.6 



X 10"^ 


U6 











The .tr^sducer systems should he phase synchronized, and 
calihratiohs should he recorded at the beginning and end of 
each day's testing. 

A sufficient humher of transducers should he available so that 
measurements can he made in and adjacent to the jet flow stream. 



h. TASK III - DATA SCALING 


U.l Scaling Procedures 

A literature search (Appendix B) failed to produce a procedure that 
appears to he more representative of current technology for scaling model data 
to full-scale aircraft sonic environments than the procedure discussed in the 
following paragraphs. 

Frequency scaling is accomplished hy considering the Strouhal number of 
the model flow field to be equal to the Strouhal number of the aircraft flow 
field Equation (ll). 


f D f D 


m 

m 

a a 

(^j 

1 

■ ("4. 

f 

m 


= model 

D 

m 


= model 

("j: 


= model 


noise data frequency 
nozzle diameter 
jet velocity 


( 11 ) 


aircraft sonic environment frequency 
aircraft engine nozzle diameter 
aircraft engine jet velocity 

. Equation (11 )■ can be solved for fg^ to obtain the aircraft sonic environ- 
ment frequency that corresponds to a designated model noise data frequency. 
Equation. (l2). 




( 12 ) 


The ratio (Djn/Dg^) is equal to the model scale. Therefore, when the model 
jet velocity is equal to the aircraft engine jet velocity (which is generally 
the case), the frequency is accomplished by a simple equation. Equation (l3) 

f = f X Model Scale (l3) 

a m 
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Amplitude scaling is easily accomplished when flow velocity and tem- 
perature of the model jet simulate the full-scale jet. Model sound pressure 
levels (SPLjh) need. only to he corrected for pressure transducer sensing ele- 
ment size (Ref. 33)' to obtain actual full-scale sound pressure levels (SPLa). 
Therefore, scaling, is accomplished hy using Equation (ll). 

SPL = SPL + 10 log-,. <|>m/<!) (Ik) 

a m 10 V /. 


The value of is determined from Figure 9. Distance from the nozzle exit 
to a location on the aircraft that corresponds to the model measurement loca- 
tion is given hy Equation (15)- 



R = distance from a point at the centerline of the aircraft engine 
nozzle exit plane to the sonic environment location of interest 

D = diameter of aircraft nozzle 
a 

D = diameter of model nozzle 
m 

R = distance from a point at the centerline of the model nozzle exit 
plane to the pressure transducer location 


When the model jet velocity and temperature do not simulate full-scale 
jet operating conditions, additional terms are required in the scaling 
equation. These addition terms are defined in Equation (l6). 


SPL 

a 


SPL^ + 10 log^o^+ K log 



+ 10 log^Q 



(16) 


K = Constant (80 for Vj <_ 6l0 m/sec and 
30 for Vj > 6l0 m/sec) 

The third and fourth terms in Equation (l6) must he scrutinized. Since 
the sonic environment on a panel .immersed in a jet flow stream is a combina- 
tion of . hydrodynamic and acoustic pressure fluctuations, a imified scaling 
equation must account for both phenomena at all points in and adjacent to the 
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jet flow field. Such an equation has not been developed to date. There- 
fore, Equation (l6) is a provisional equation that has been defined in order 
to scale the SCAT 15-F data to full-scale supersonic transport engine oper- 
ating conditions. The velocity and density terms are based on methods for 
predicting acoustic power of free jets as defined in References 35 and 36. 

The basic assumption made to derive Equation (l6) was that acoustic pressure 
fluctuations have a greater impact than hydrodynamic pressure fluctuations on 
the sonic environment of a panel immersed in the flow field of a high temper- 
ature supersonic jet. If the hydrodynamic flow field is considered to have a 
predominant impact on the sonic environment, the sound pressure levels will 
be a function of dynamic pressure (q). The relation normally used is 
SPL=20 log q + constant. Therefore, since q=l/2 pV?, the scaling equatdon 
will be in terms of velocity to the fourth power. Reference 37 shows that 
the sonic environment of surfaces immersed in a jet flow stream increases 
according to V^. The relative importance of acoustic and i hydrodynamic pressure 
fluctuations will not be pursued fiirther in this report. However, a computer 
program has been developed (Appendix C) based on Equation (l6) and is 
subsequently evaluated in Section U.2. 


^ . 2 Comparison of Scaled Model Data and Full Scale Aircraft Date 

The literature search (Appendix B) did not result in finding a corre- 
sponding set of model and full-scale aircraft data that is pertinent to sonic 
environment in a supersonic jet flow stream. However, subsonic jet data was 
available for the following three aircraft. 

• S-3A (Refs. 38 and 39) 

• L-1011 (Refs. 39 and itO) 

• V/STOL (Ref. itl) 

The full-scale sonic environment data for the V/STOL aircraft were 
obtained from unpublished data that were provided by Langley Research ■ Center . 

To compile the model test data in a format for input to the computer 
program, the- following information must be available. 

• Model Data: scafac, ttmr, ptmi, dmi , inixm, rmi , trim fm, sp£m 

• Aircraft Data: ttar, pra, wa, mixa, tria 
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data 


.. 2 


• Atmospheric Data; psi 

scafac = model scale factor 

ttmr = total temperature of model jet - deg R 

ptmi = total pressure upstream of the nozzle Ih/in'" 

dmi = model nozzle diameter (exit) - in. 

mixm = fuel/air mixture of model jet 

rmi = distance from model nozzle exit to pressure transducer - in. 

trim =■ diameter of pressure transducer sensing element (model test) - in. 

fm = model frequency band center frequency - Hz 

sp5,m = sound pressure levels corresponding to fm - dB 

ttar = total temperature of engine jet - deg R ■, 

pra = engine nozzle pressure ratio 

wa = engine exhaust flow'- Ib/sec 

mixa = engine exhaust fuel/air ratio 

tria = diameter of pressure transducer sensing element (Full scale 
test) - in. 

2 

psi = atmospheric pressure - Ib/in 

Several problems were found to exist in compiling the required input 
These included. 

• Sufficient information was not reported with regard to propulsion, 
geometry, and instrumentation. 

• Decisions had to be made regarding the relative importance of the 
core and bypass jet streams 

• Spectra that was expressed in terms of spectral density had to be 
converted to third octave band data. 

• Range of the frequency spectr\un was too small 

• Validity of test data could not be assessed. 
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The data compiled in Table T were determined to be the best available 
and were used for scaling the model data to full-scale conditions. 

Figure 21' shows location of the pressure transducers for the S-3A. The 
locations were similar for the model and aircraft. Figures 22, 23 and 2k 
show the measured model data, measured aircraft data and scaled sonic envi- 
ronment for pressure transducer locations 1 , 2 and 3, respectively.. 

Figure '25 ' shows the transducer location for the L-1011 tests and Fig- 
ure 26 shows the comparison between scaled and measured noise. 

Figures 27 and 28 give the geometrical and noise data for a V/STOL air- 
craft with lower surface blowing of the flaps by the jet stream. 

The S-3A and L-1011 model jets sim\ilated fxill-scale aircraft engine 
characteristics. Therefore, velocity and density terms in Eq.uation (I 6 ) 
have a minimal effect on noise scaling. Consequently, scaling of the model 
data is accomplished according to Equations (13) and (l^) for the frequency 
and amplitude respectively. Figures 22 through 2k and Figure 26 show that 
the scaled noise spectrum has approximately the same shape and magnitude as 
the' measured model spectrum and has a frequency shift that is proportional 
to the model scale factor. Pressure transducers 1 and 3 on the S-3A air- 
craft (Figure 21) and the transducer on the L-1011 aircraft (Figure 25) are 
believed to be in the jet flow stream whereas the number 2 S-3A transducer 
is outside of' and adjacent to the flow stream. Observation shows that the 
scaled model spectrum for the No. 2 S-3A transducer agrees reasonably well 
in shape and magnitude with the measured full-scale data whereas the scaled- 
model spectra in Figures 22, 2k and 26 are not representative of the full- 
scale measured spectra. Comparison of the scaled spectra with full-scale 
measured data does not show any specific trend regarding shape, peak fre- 
quencies or magnitude. 

Scaling of the V/STOL transducer data (Figure 28) for the transducer 
location shown in Figure 27 involved accounting for velocity and density 
differences in model and full-scale jets. Figure 28 shows that the scaled 
model data is significantly higher than the full-scale measured data. 

Discrepancies between the scaled model data and full-scale measured data 
may be attributed to several factors. These factors include: 

• The models and full-size engines had coaxial nozzles and bypass jet 
parameters were used to scale the data 

• Noise measurements were made by different investigators which used 
different types of instriamentation and had different test 
environments . 

• The compatible model and full-size aircraft data that were available 
for evaluating the prediction method were for subsonic jets whereas 
the scaling equation was developed for supersonic jets. 


52 



TABLE 7- COMPUTER PROGRAM INPUT DATA 


S-3A 

L-1011 

V/STOL 

LOC 1 

LOC 2 

LOC 3 




1 

0.143 


0.05 

0.185 



537 



537 


537 




20.8 


22.2 

22.2 



4.2 


2!88 

6.95 



0 


C 


0 


8 


18 


30 

11.5 

30.5 



0.25 


0.25 

0.218 



593 



620 


1,515 




1.42 


1.49 

1.5 



411 



1,159.5 

333.5 



0 


C 


0.033 



0.25 


0.25 

0.5 



14.7 


14.7 

14.7 

fm 

splm 

fm 

splm 

fm 

splm 

fm 

splm 

fm 

splm 

250 

142.7 


132.2 


130.7 

1,250 

126.4 

315 

144.1 

315 

143.6 


132.3 


127.5 

1,600 

126.2 

400 

144.9 

400 

145.1 


132.5 


125.4 

2,000 

125.3 

500 

145.6 

500 

145.6 


130.2 


121.8 

2,500 

125.3 

630 

145.4 

630 

146.4 


128.2 


120.4 

3,150 

124.2 

800 

145.1 

800 

146.3 


126.2 


119.5 

4,000 

123.2 

1,000 

144.9 

1,000 

146.8 


125.0 


120.0 

5,000 

122.0 

1,250 

144.7 

1,250 

147.1 


124.4 


120.5 

6,300 

120.2 

1,600 

143.8 

1,600 

146.2 


123.3 


120.4 

8,000 

119.0 

2,000 

143.9 

2,000 

145.5 


122.1 


120.0 

10,000 

117.0 

2,500 

. 143.0 

2,500 

144.5 


120.0 


119.5 

0 

0 

3,150 

142.3 

3,150 

142.8 


118.2 


117.9 

0 

0 

4,000 

141.0 

4,000 

141.5 


117.0 


117.8 

0 

0 

5,000 

139.0 

5,000 

139.8 


116.2 


116.4 

0 

0 

6,300 

137.8 

6,300 

138.6 


115.0 


115.3 

0 

0 

8,000 

136.2 

8,000 

136.6 


114.0 


113.4 

0 

0 

10,000 

133.0 

10,000 

134.2 


112.4 


110.5 

0 

0 

12,500 

129.2 

0 

0 


0 


0 

0 

0 

16,000 

125.5 

0 

0 


0 


0 

0 

0 

20,000 

121.6 

0 

0 


0 


0 

0 

0 

25,000 

116.3 

0 

0 


0 


0 

0 

0 

31,500 

112.0 

0 

0 


0 


0 

0 

0 

40,000 

110.0 






































THIRD-OCTAVE BAND LEVEL IN DB RE 2 x 10‘° N/m 


THE COMPUTER PRINTOUT SHOWING THE MODEL SPECTRA, 
FULL SCALE SPECTRA, AND SCALED SPECTRA FOR FIGURES 
22, 23, 24, 26 AND 28 ARE GIVEN IN APPENDIX C. 



Figure 22. S-3A Spectra - Location 1 
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THIRD-OCTAVE BAND LEVEL IN DB RE 2 x lO'^ N/m^ 
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TEST CONDITIONS: TABLE 7 
COMPUTER OUTPUT: TABLE 10 
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Figure 24. S-3A Spectra - Location 3 
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RB-211 

ENGINE 


Figure 25. L-1011 Measurement Location 







THIRD-OCTAVE BAND LEVEL IN DB RE 2 x 10*° N/m 



Figure 26. L-1011 Spectra 
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ENGINE 


Figiire 27. 


V/STOL Meastirement Locations 








^.3 SCAT 15 -F Scaled Model Data 

Inability to accurately scale the data in Section h.2 may be due to the 
complex confluent flow field ^ of the coaxial jets. Therefore, even with the 
limitations of the SCAT 15-F model test data (see Section 3.3), scaling of 
the data to full size aircraft dimensions should provide some insight as to 
the severity of the sonic environment for a SST with over-the-wing engine 
installations. Figure 29 gives a comparison of SCAT 15-F measured model data 
with scaled spectra for a realistic SST engine. The spectra plotted are for 
run nianber 31P transducer 10. Model data and scaled spectra are given in 
Table 15 of Appendix D along with model and full scale engine operating char- 
acteristics. Table I 5 shows model velocity and full-scale engine velocity 
to be hkO m/s ( 1 U 5 O .5 ft/s) and 845 m/s (2766 ft/s), respectively. This 
differential velocity has a significant impact on the scaled noise level, see 
Equation ( 16 ). The computed value, because of decreased density for the 
full-scale engine jet, is shown in Table 15 to be -6.68 dB. Comparison for 
other run numbers and transducer locations could be made by plotting values 
given in Tables 16 through 30. 

Figure 30 shows scaled OASPL contours on the wing surface of an arrow 
wing for Mach 1.5 nozzle cold jet test conditions. This is the same contour 
plot that is given in Figure 16 . Maximum OASPL values occur approximately 
22 meters downstream of the nozzle exit. Since the scaled-nozzle diameter is 
1.69 meters, the ratio of downstream distance to nozzle ' diameter is 13. 

OASPL contour plots were not made, for the scaled SST engine operating 
conditions. However, Figure 31 gives a point-by-point comparison of scaled 
Mach- 1.5 nozzle data and realistic SST engine OASPL values-. Differences in 
the Mach 1.5 scaled OASPL values and realistic SST engine OASPL values range 
from 6 to 20 dB because of sensitivity of the scaling method to spectrum 
shape. The peak sonic environment level for the wing structure is shown to. 
have an OASPL value of 18 I dB. 
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THIRD-OCTAVE BAND LEVEL IN DB RE 2 X 10 ° N/m 
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Figure 29- Comparison of SCAT 15-F Measured Model 
Data with Scaled Spectra for a 
Realistic SST Engine 






Figure 31 . Point t>y Point Comparison of Scaled Mach 1.5 Nozzle Dat 
and Realistic SST Engine OASPL Values 



5 . CONCLUSIONS 


Accuracy of the prediction method used for scaling the SCAT 15-F test 
data to full-size supersonic transport dimensions has not been established. 
However, scaling of the SCAT 15-F data by considering the bypass flow condi- 
tions of an SST duct-burning turbofan engine concept,, results in noise levels 
of up to 180 dB on. a surface immersed in the jet flow Stream. It is not 
inconceivable to believe that noise levels of this magnitude are possible for 
jet velocities that are on the order of 1000 m/s as are typical of duct- 
burning turbofan engines. 

Sufficient test data pertaining to the sonic environment of a surface 
immersed in a supersonic jet flow stream does not exist to compile a data 
bank for determining validity of model scaling procedures. However, scaling 
of the S-3A, L-1011 and a V/STOL subsonic jet model data to full-size air- 
craft conditions and making comparisons with measured aircraft data show 
discrepancies of up to approximately 20 dB in some of the spectra one-third- 
octave bands. These results indicate that additional studies are needed to 
evaluate scaling methods and model test procediares. 

The relation between sonic environments on a structure that are associ- 
ated with hot and cold jet flow fields has not been determined. Hot jet flow 
sonic environment measurements have not been obtained because pressure trans- 
ducers have not been developed that are known to provide accurate measurements 
in a hot jet flow stream. Therefore, significant improvements in. high- 
temperature pressure transducers are needed. 

Finite size of a transducer sensing element limits its space resolution 
of a pressure field. The error becomes progressively larger as the value of 
the quantity coR/uq increases . Therefore, values of the circular frequency 
should be as low as possible. This can be accomplished by increasing the 
model size so that the model frequency range required is made smaller. Also, 
the sensing element of the transducer should be as small as possible. The 
sensitivity of the transducer decreases as the sensing element becomes 
smaller. Consequently, a compromise must be made between sensitivity verses 
spatial resolution characteristics. 

If na,rrow-band cross-correlation functions are required for making modal 
analyses, the transducer systems must be phase synchronized. Otherwise, an 
initial phase shift will result in a time shift of the entire function. 

Scaling of data from coaxial jets is more complex than scaling the data 
for simple circular nozzles. Therefore, it appears that simple circular 
nozzle tests will provide a better basis than coaxial nozzle tests for inves- 
tigation and refining of supersonic jet noise scaling procedures. 
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APPENDIX A 


DERIVATION OF RESPONSE SPECTRAL DENSITY EQUATION 


Determining the response of a structure that is subjected to a random 
excitation force consists of three major steps. They are: 

• Expression of the random excitation force in terms of the excitation 
spectral density 

• Determining the receptance of the system 

• Expressing the response spectral density in terms of receptance 
and excitation spectral density 

The computation procedure will first be developed for a simple spring mass 
system and then extended to the general procedure for a structure. 


A.l Simple Spring-Mass System 

Figure 32 shows the steps required for determining the response of a 
simple spring-mass system. Plot a is a time history of the excitation force. 
Plot b is a typical excitation spectral density curve, Plot c expresses the 
receptance in terms of the absolute value of the receptance squared and Plot d 
shows the response spectral density. Computation methods for obtaining plots 
b, c and d are given below. 

Excitation Spectral Density . - The excitation spectral density can be 
obtained by relating the spectral density to the autocorrelation function by 
use of the Fourier Transform, Equation (IT)- 



Rp(x)E 


-io)T 


dT 


(IT) 


Sp(m) = Excitation spectral density 
Rp(x) = Autocorrelation function 
0 ) = Circular frequency 

T = Delay, time 
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P(t) 




(a) 

EXCITATION 
TIME HISTORY 



(b) 

EXCITATION 
SPECTRAL DENSITY 



(c) 

RECEPTANCE 
OF STRUCTURE 



(d). 

STRUCTURAL 
RESPONSE 
SPECTRAL DENSITY 
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Figure 32. Steps Involved in Determining the Structural 
Response Caused hy a Random Excitation Force 



The autocorrelation function can be determined from the random excitation 
signal by using an autocorrelation function analyzer to determine the time 
average of the product of the signal at times t and t + t. Equation (l8) 

R(t) = <P(t) P(t + t)> (i8) 

X ; = Delay time 

. P(t) = Value of signal at time t 

P(t + x) = Value of signal at time (t + x ) 

Receptance of the System. - The equation of motion for a simple spring- 
mass system (Figure 33) is given by Equation (19)- 


mx+cx+kx=P(t) (19) 

The response to a transient force can be determined from the receptance 
by using a Fourier integral technique, but it is usually more convenient to 
make use of a convolution integral which expresses the response in terms of 
the response to a unit impulse. An impulsive loading can, be expressed mathe- 
matically by Equation (20). 

P(t) = I6(t) (20) 

P(t) = exciting force 
I = magnitude of the impulse 
(S(t) = dirac delta function 

The response to the impulse loading is given by Equation (21). 

x(t) - W(t) I (21) 

x(t) = System response 

W(t) = Response to a unit impulse 

The solution of the equation of motion when P(t) = 16 (t) is given by 
Equation (22). 


X 


I -(c/2m)t 
— — ' e ' 
mu)^ 
a 


sin to.t 
a 


(22) 
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P(t) = EXCITING FORCE 
x(t) = RESPONSE OF SYSTEM 
m = ELEMENT OF MASS 
c = DAMPING COEFFICIENT 
k = SPRING CONSTANT 


Figure 33. Simple Spring Mass System 



t 


Figure 3U- Continuous Pressure Loading Plot 
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where 


OJ 


d 



2 


= Damped natural frequency 


Substitution of Equation (22) into Equation (21) and solving for W(t) 
results in Equation (23). 


W(t) 


1 -(c/2m)t 

— ^ — e 

mil). 


sin 


0),t 

d 


(23) 


If the loading is continuous instead of a single impulse (see Figure 3;4) , 
the area under the loading-time curve can he divided into impulsive loadings 
P(t)6t and the response at any time t is given hy the convolution integral, 
Equation (24). 


x(t) 


I 


W(t - T ) P( T ) dT 


(24) 


The value of the convolution integral is not changed hy a shift in time. 
Therefore, the response can also he determined hy Equation (25). 


■I 


CO 

x(t) - I W(t) P(t - t) dT 
o 


(25) 


If P(t) = P^ e^^^ the response is given hy Equation (26) 


x(t) = 


P e' 
o 


i(jjt 


W(T)ie dx = 6(ia)) 


(26) 
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where 


6(103) = 


W ( T ) e dT 


6(io3) = Receptance of the system 


Equation (26) will he used subsequently for determining the response 
spectral density. The value of the receptance can be determined by determin- 
ing the particular integral of Equation (27)- 


mic + cx + kx 


= P 

o 


(27) 


Divide Equation (27) by m then let c/m = 2^0) and k/m = 
Equation (28). 


03 to obtain 
n 


2 Po lojt 

X + 2Co3 X + 03 X ,e 

n ■ n m 


(28) 


The solution of Equation (28) is x = x^e^*^ . Therefore, by taking the 
first and second time derivatives of the solution and substituting the 
results into Equation (28), Equation (29) is derived. 


2 2 
(-03 + i2?0) 03 + 03 ) 

n n 


X 

o 


p 

o 

ra 


(29) 


Solving for x^ gives Equation (30). 


X 

o 


P 

o 

2 2 

m(o) - 03 + iC03 03 ) 

n n 


(30) 


The required response. Equation (31), is determined by substituting the 

value for x into the solution, 
o 
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1 ^ - io)t 

2 2 

m(o) - 0 ) + i^(jj (u) 

n n -J 


The term in the brackets is the reeeptance of the system. Equation (32). 


a(itio) = 


2 2 

m(w - 0 ) + i6(jj w) 

n n 


The ordinate of the plot in- Figure 32-c'is the square of the absolute 
valve of the reeeptance. 

Response Spectral Density. - The response autocorrelation function is 
given by Equation (3S)- ^ 


R^(t) = <x(t) x(t■^T)> 


Substitution x(t') from Equation (25) gives; 


x(t) = I W(t^) P(t-T^) dx 


x(t+x) = W(t^) P (t-^T-T 2 ) dx 


Substitution of Equation (3^) into .Equation (33) gives; 


R (t) =1 W(x ) I W(x ) <P(t-x );'P(t-Hx-X 2 )>dX 2 dx, 
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By changing the origin of t: 


J fOO CO 

I W(t^) 

r, 


R^(t) = I W(t^) I <P(t) P (t+T^-T2T)>dT2<iT^ 


(36) 


The time average <P(t) P (t + + t)> is the autocorrelation func- 

tion of the exciting force. Therefore: 


oo oo 


'• r 

= W(x^) 


B^(t) = j W(x^) I ^{ t ^) Rp(t^ ~ '^2 


(37) 


The relationship "between the response and excitation spectral densities 
can he obtained from Equation (37) hy making use of the Fourier transform rela- 
tionship between the autocorrelation function and spectral density, Equa- 
tion ( 38 ). 


s^{») 



(t)e 


-itOT 


dx 


(38) 


Substitution of Equation (21) into Equation (38) gives: 


-1U)X 


r°° r°° 

SM = / [j W(t^) I W(t^ Rp (Tp-Tp+T) dTpdrJ e- 

-OO •'o *^0 


(39) 


■r 


W(x^)e^“’^l dT 




fvu- 

•' — on 


i2+T)e dCx^-x^+T) 
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The separate factors of the above expression are a*(iw), a(.i(o) and 
Sp(a)) . a*(io)) is the complex conjugate of the receptance and Sp(o3) is the 
excitation spectral density. Therefore, the response spectral density is 
given by Equation (UO). 


S^(o)) = a*(ico) a (iw) Sp(ai) 



(J+0) 


This is the plot shown in Figure 31-d. 


A. 2 Response of a Structure 

A detailed development of the generalized structural response spectral 
density Equation is given in References and ^43. Therefore, only the basic 
equations will be given herein. The response spectral density of any point 
(x-^) will be detennined for distributed pressures at points Xj^ and Xg for a 
beam of length £. This development - shows all the essential features for more 
complicated systems. The beam problem requires that the displacement at any 
point be expressed in terms of the normal modes and normal coordinates as 
given in Equation ,(i+l ) . 


w(x,t) = Z v^(x) C^(t) (i+l) 

r 

w(x,t) =. Displacement at any point 

w (x) = Normal mode of the beam 

r 

5^(t) = Normal coordinate 

It should be noted that the normal coordinates are obtained by transform- 
ing the generalized coordinates of a system by the normal mode matrix. When 
the equations of motion are expressed in terms of the normal coordinates the 
system is inertially and elastically uncoupled. Therefore, the equations do 
not have to be solved simultaneously. 

For a body ^f any shape the position at anyj.point is described by a 
position vector P and the displacement vector w(P,t) is given by Equation {k2) . 


w(P,t) = Z w (P) C (t) 
r r 
r - 


(U2) 
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Therefore, the only difference between the beam equations given below 
and the equations fo£ a body of any shape is that the coordinate x is used 
instead of a vector P. The steps involved in determining the generalized 
response spectral density equation are the same as for the simple spring-mass 
system. However, the crosfs correlation functions and cross spectral densities 
must now be determined. Also, the receptances of the system must be defined 
in terms of the normal modes. 


Excitation Spectral Density . - The excitation cross spectral density is 
obtained by relating it to the cross-correlation function by use of the Fourier 
Transform Equation (i^3). 


Sp(x^, Xp; m) 



(3h) 


The cross-correlation function is: 


Rp(xA, Xp; t; 


t+T)dx3> 


(i+M 


Sp(xA, Xp-, m) 


Rp(xA, Xp, t) 

P(x^,t) 

P(Xp, t+r) 


= excitation cross spectral density at any two 
points x^ and Xp 

= cross-correlation function for points x and x.^ 

A 

= distributed pressure at point x^ at time t 
= distributed pressiire at point x_ at time (t + t) 

ID 


Receptance of the system. - The receptance of the system at point x^ for 
a load at x^ in terms of the normal modes is given by Equation (^+5)- 


a = Z 


2 2 

M (o) - CO + i^o) co) 

r r r 




The receptance at x^ for a load at Xp is; 


a = E 

X., ,x„ s 


w (x ) w (x ) 
si s B 


l’'^B “ M (co^ - cj^ + iCco co) 

s s s 
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W (x ) , V (x ) 
r 1 s 1 






•= modal deflections at x^ 
= modal deflection at 
= modal deflection at x^ 


M 

r 



w^(x) mdx 


= generalized mass for mode r 


M = I 

' J 


w (x) mdx = generalized mass for mode s 
s 


0 ) , m = modal frequencies 

r 'S 


? = viscous damping factor 

Response Spectral Density . - The spectral density at x^ can now be 
expressed in terms of the receptances. and excitation cross-spectral density 
by Equation (46). 


S (x m) 
w 1 r s 


Z I 

u 


w^(Xi)w^(xi) 4 I w^(x^)w^ (xg)Sp(x^,x^;u,) dx^dx^ 

2 2 2 2 " 
M (oj - 0) + iCo) w) M (o) - 0) + iC(jJ w) 

r r r X s s 


(46) 


This equation is easily changed to Equation 6 in Section 2.3 by: 

• Replacing x, with x,y 

• Changing integration limits to A 

• Letting x^ = ^ 

• Letting x.^ = p 

r> 

• Replacing dx. and dx^ with dA 

• Representing modal values w^ and w^ by f^ and f^ respectively 

• Changing symbol for damping factor from 5 to 6 to 6 ^ and 6^ 

respectively ^ 
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APPENDIX B 


LITERATUEE SEARCH 


A search was made of the following files via DIALOG - the Lockheed- 
California Company on-line information retrieval system - to locate informa- 
tion relevant to " Acoustic Loads on a Panel immersed in a Jet Flow Stream. " 

e NTIS - File of Government Reports: The file contains over U00,000 

abstracts of research reports from over 2^0 government agencies, 
including NASA. The file dates from 196^+. 

• ENGINEERING INDEX - Publications of engineering organizations: The 

file contains approximately 360,000 citations and abstracts from 
3,500 Journals. The file dates from 19T0. 

• ISMEC - Mechanical engineering and engineering management data base: 
The file includes about 30,000 items and dates from 1973. 

A search was made by the NASA Scientific and Technical Information Facility. 

The search also included: 

• Journal of the Acoustic Society of America Index - 1971 to present 

• Lockheed California Company Central Library Catalog 

o , Applied Mechanics Services 197^ to Sept 1975 

• Applied Science and Technology Index - 1971 to Sept 1975 

• AGARD Index, - 1971 thru 1973 

• Shock and Vibration Digest - 1971 to Sept 1975 

Test Data directly related to "Acoustic Loads on a Panel Immersed in a Jet 
Flow Stream" were not located. However, many publications were located which 
provide pertinent information for development of sonic environment analysis 
methods (see report references). 
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APPENDIX C 


SCALING COMPUTER PROGRAM 


A listing of the noise scaling computer program is given in the following 
pages. The program is written in terms of the International Business 
Machine (IBM) conversational program system (CPS) CPS PL/l language. The. 

CPS PL/l language can he regarded as a modified subset of the full set PL/l 
language . 


Input to the program is: 

I scafac I ttmr | ptmi | mixm | dmi | rmi | trim 

I fm I splm I n I 

I ttar I pra ( wa j mixa | triaj 

I Psi I ■ 1- 


These input symbols are defined in the program symbol list that follows 
the program listing. Input req,uired is in terms of the English ' system of 
units . 
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noise' scaling computer program 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 

10 . 

11 . 

12 . 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20 . 

21 . 

22 . 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 


42. 

43. 


DECLARE fm(24) DEC ( 6), fa (24) DEC ( 6) ,asp1 a ( 24) DEC(6), 
splm(24) DEC(6); 

DECLARE tabsc(12) DEC ( 6) , tord( 12 ) DEC ( 6 ) ,mspl a ( 24 ) DE 
C(6); 

DECLARE scm(24) DEC ( 6) , tc( 24) DEC ( 6 ) , sea ( 24) DEC(6); 
DECLARE gam ENTRY EXT KEY(wag); 

GET LI ST( scafac, ttmf/ ptml ,ml xm, dml / rml , t r Im); 

GET L I ST( fm, spl my n ); 

GETLIST(ttaryprayWaymlxaytrfa); 

GET L ISTCpsI ); 

GET LIST(start); 
dmf=dmt / 12; 
daf=dmf/scafac; 
rmf=rml/12; 
raf=rmf/scafac; 
prm“ptml /ps t ; 
cvycvmycva= .98; 
aml=.785*dmt**2; 
amf=aml/144; 

CALL gamCml xa, ttar, praygamaa, tj ar); 

CALL gam(mlxmy ttftiry prmygamamy t j mr); 

vm=sqrt (64 . 4*( gamam/ (gamam-1 ))*53.3*ttmr*( 1-1/prm** ( ( 
gamam-1 )/gamam) ) ); 

va»sqrt(64 . 4* (gamaa/ (gamaa-1 ))*53.3*ttar*(l-l/pra**(( 
gamaa-l)/gamaa) ) ); 

rhom^ps I *144/(53 . 3* ttmr)*prm**( (gamam-1 )/gamam); 

rhoa=psI *144/(53.3*ttar)*pra**( (gamaa-1) /gamaa); 

wm=rhom*amf*vm*cv; 

aa=wa/ (cv*va* rhoa ) ; 

fnm=cv* (wm/32 .2 )*vm; 

fna=cv* (wa/32 . 2 ) *va; 

sd=10*l ogl0( rhoa/ r horn); 

127: prlnt“l; 

ucm“ . 62*vm; 

DO 1=1 TO n; 

omegam=6.28*fm( I ); 
abscm=omegam* ( t r Im/24)/ ucm; 

DO j=l TO 11; 

IF abscm>=tabsc(j )&abscm<"tabsc( j +1) THEN GO TO 1 
15; 

END ; 

115: phi rm= (abscm-tabsc( j ))/(tabsc(j+l)+tabsc(J))*(tord( 

J +1) -tord( j ) ) +tord(J ); 
scm( 1)=-I0*logl0(phlrm); 
tc( 1 )»scm( I ) +sd; 
fa( 1 )»fm( I )*(va/vm*(dmf/daf ) ); 

IF va<=2000 THEN aspla ( 1 )=spl m( I ) +scm( i ) +sd+80* logl 
0(va/vm); ELSE aspl a ( I )=spl m( 1 ) +scm( I ) +sd+165 . 051+3 

0* 1 oglO ( va ) -80*1 oglO(vm); 

END ; 

uca= .62*va; 
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44. 

45. 

46. 

47. 

48. 

49. 

50. 

51. 

52. 

53. 

54. 

55. 

56. 

57. 

58. 

59. 

60. 
61. 

62. 

63. 

64. 

65. 

66 . 

67. 

68 . 

69. 

70. 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 

79. 


NOISE SCALING COMPUTER PROGRAM - Continued 


140: 


DO !■! TO n; 

omegaa»6 . 28*fa ( I ); 
absca=omegaa*( tr la/24)/uca; 

DO j«l TO 11; 

IF absca>“tabsc( j )&absca<“tabsc( j ♦!) THEN GO TO 1 
40; 

END ; 

phi ra=(absca-tabsc(j >)/(tabsc(j +l)+tabsc(j ))*(tord( 
j+l)-tord(j ))+tord(j ); 

. sca( I )=10*1ogl0(ph?ra); 

IF l<«n THEN msp1a( I )=aspla( I )+sca( I ); ELSE mspla=0 

• 

END ; 

PUT IMAGE(pr!nt)(110); 

110: IMAGE; 

SONIC ENVIRONMENT SCALING PROGRAM 
PUT LIST(lf(2)); 

PUT IMAGE(prInt)(117); 

117: IMAGE; 

MODEL AIRCRAFT 

PUT LISTC '); 

PUT IMAGE(ttmr,ttar)(118); 


118: IMAGE; 

tt 

PUT IMAGE(tjmr,tjar)(119); 
119: IMAGE; 


ts 


1 25 : 


pr 


153 


w 


121 : 


155 


ml X 


cv 


PUT IMAGE(prm,pra)(125); 
IMAGE; 

PUT IMAGE(wm,wa)(133); 
IMAGE; 

PUT I MAGE (ml xm,mt xa ) ( 1 21 ) ; 
IMAGE; 


PUT I MAGE (cvm,cva ) ( 1 35 ); 
IMAGE; 


PUT IMAGE(dml,daf)(123); 

123: IMAGE; 

d I a . — 1 n 

PUT IMAGE(aml,aa)(130); 

130: IMAGE; 

area . — sq In 

PUT IMAGE(rhom,rhoa)(132); 
132: IMAGE; 

rho — . — 

PUT IMAGE(gamam,gamaa)( 1 26); 
126: IMAGE; 

gama — — 


deg R 
deg R 


1 b/sec 
#f/#a 


ft 

sq ft 
Ib/sq ft 
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NOISE SCALING COMPUTER PROGRAM - Concluded 


80. PUT IMAGE(vm,va) (1 22); 

81. 122: (MAGE; 

vel .- .- ft/sec 

82. PUT IMAGE(ucm,uca)(136); 

83. 136: IMAGE; 

conv ft/sec 

84. PUT IMAGE(rml\,raf)(124); 

85. 124: IMAGE; 

dlst . -- In . — ft 

86. PUT IMAGE(fnm, fna)(134); 

87. 134: IMAGE; 

Fn .- .- lbs 

88. PUT IMAGE(trlm,trla)(138); 

89. 138: IMAGE; 

trdia “. In 


90. PUT LIST(lf(2)); 

91. PUT IMAGE(sd)( 137); 

92. 137: IMAGE; 

Density Correction = . — dB 

PUT LIST(lf(2)); 

PUT IMAGE(prlnt)(113); 

IMAGE; 

splm scm tc fa 

sea mspla ; 

PUT IMAGE(prlnt)(139); 

IMAGE; 

(dB) (dB) (dB) (Hz) 

(dB) (dB) 

PUT LIST( • '); 

DO 1=1 TO n; 

PUT IMAGE(fm( I ),splm( I ),scm( I ) , tc( I ) , fa ( I )/ aspl a ( I ) 
, sea ( 1 ) ,mspl a( I ) ) ( 1 14); 

IMAGE; 


102. 


END ; 

103. 


S“0; 

104. 


DO 1=1 TO n; 

105. 


s=10**(mspla(l)/10)+s 

106. 


END ; 

107. 


oaspi =10*1 ogl0( s) ; 

108. 


PUT LISTdfd)); 

109. 


PUT IMAGE(oaspl )(120); 

110. 

120: 

IMAGE; 


OASPL= 

.- dB 

111. 


PUT LISTdfdS)); 


93. 

94. 

95. 

aspi a 

96. 

97. 

(dB) 

98. 

99. 

100 . 


113: 

fm 


139: 

(Hz) 


101 , 


114; 
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COMPUTER PROGRAM VARIABLES 

Units 

aa 

= 

Aircraft engine fully expanded jet area 

ft^ 

absca 

= 

Abscissa of finite size correction curve 
(Figure 9) for aircraft noise 

— 

abscm 

= 

Abscissa of finite size correction curve 
(Figure 9) for model noise 

— 

ami 

= 

Area of model jet 

in2 ■ ■ 

amf 

= 

Area of model jet 

ft^ 

aspla 

- 

Aircraft noise overall soimd pressure level 

dB 

cva 

- 

Nozzle velocity coefficient of aircraft engine 
nozzle 

— 

cvm 

= 

Nozzle velocity coefficient of model engine 
nozzle 

ft/sec 

daf 

= 

Diameter of aircraft engine nozzle exit 

ft 

dmf 

= 

■ Diameter of model nozzle exit 

ft 

dml 

= 

Diameter of model nozzle exit 

in. 

fa 

= 

Aircraft noise frequency (one-third octave) 

Hz 

fm 

= 

Model noise frequency (one-third octave) 

Hz 

fna 

= 

Aircraft engine jet thrust 

lb.:- 

fnm 

= 

Model jet thrust 

lb- 

gam 

- 

Subroutine for computing values of ratio of 
specific heat 

— 

gamaa 

— 

Ratio of specific heat for aircraft engine fully 
expanded jet 


gamam 

--- 

Ratio of specific heat for model engine fully 
expanded jet ' 

— 

i 

= 

Do loop counter 

— 
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.COMPUTER PROGRAM VARIABLES - Continued 





Units 

j 

= 

Do loop counter 

— 

mixa 

= 

Fuel/air mixture of aircraft jet 

Ib-fuel/lb-air 

mixm 

= 

Fuel/air mixture of model jet 

Ib-fuel/lb-air 

mspla 

= 

Model noise overall soimd pressure level 

dB 

n 

= 

Humber of one-third octave bands 


omegaa 

= 

Aircraft noise one-third octave band circular 
frequency 

rad/ sec 

omegam 

= 

Model noise one-third octave band circular 
noise 

rad/sec 

phira 

= 

Transducer finite size correction for aircraft 
noise 

dB 

phirm 

= 

Transducer finite size correction for model noise 

dB 

pra 

= 

Aircraft engine nozzle pressure ratio 

— 

prm 

= 

Model nozzle pressure ratio 

— 

print 

= 

Dummy variable used to facilitate printout 

— 

psi 

= 

Atmospheric pressure 

Ib/in^ 

pta 

, - 

Total pressure of aircraft jet 

Ib/ft^ 

ptmi 

= 

Total pressure of model jet 

Ib/in^ 

{ 

raf 

- 

Distance from aircraft engine nozzle exit to 
transducer 

ft 

rmf 

= 

Distance from model nozzle exit to transducer 

ft 

rhoa 

= 

Fully expanded aircraft jet density 

slug/ ft^ 

rhom 

= 

Fully expanded model jet density 

s lug/ft ^ 

rmi 

= 

Distance from model nozzle exit to transducer 

in. 

sea 

= 

Transducer finite size vorrection for aircraft 
noise 

dB 



COMPUTER PROGRAM VARIABLES - Concluded 





Units 

scafac 

= 

Model scale factor 

dB 

scm 

= 

Transducer finite size corrections for model noise 

dB 

sd 

= 

Density correction 

dB 

spla 

= 

Aircraft noise one-third octave band sound 
pressure levels 

dB 

splm 

= 

Model noise one-third octave band sound 
pressure levels 

dB ■ 

start 

- 

Dummy variable used to initiate start of analysis 
and output 

. ■ — 

tabsc 

= 

Abscissa values of Figure 9 

— 

tc 

= 

S\un of finite size and density corrections 

dB 

tjar 

= 

Static fully expanded aircraft Jet temperature 

deg R 

tjmr 

= 

Static fully expanded model Jet temperature 

deg R 

tord 

= 

Ordinate values of Figure 9 


tria 

= 

Diameter of transducer sensing element (aircraft 
test ) 

in. 

trim 

= 

Diameter of transducer sensing element (model 
test) 

in. 

ttar 

• = . 

Total temperature of aircraft engine Jet 

deg R 

ttmr 

= 

Total temperature of- model engine jet 

deg R 

uca 

= 

Convection velocity of aircraft engine Jet 

ft /sec 

ucm 


Convection velocity of model Jet 

ft/ sec 

va 

= 

Aircraft engine fully expanded Jet velocity 

ft/ sec 

vm 

= 

Model fully expanded Jet velocity 

ft/ sec 

wa 

= 

Aircraft engine exhaust flow 

lb/ sec 

wm 

- 

Model engine exhaust flow 

lb/ sec 



APPENDIX D 


COMPUTER PROGRAM OUTPUT 

included in this appendix are the output data for: 

• S-3A (Tables 8 through 10) 

• L-1011 (Table 11 ) 

• V/STOL (Table 12) 

• -SCAT 15-F (Tables 13 through 30 ) 

The symbols used on the printout sheet correspond to those given in 
Appendix C. However for convenience the symbols for each column are defined 
below. 


1. 

fm I 

■ = 

Model frequency 

2. 

splm 

= 

Model noise one-third octave spl 

3. 

scm 

= 

Transducer finite size correction for model noise 

It. 

tc 

= 

Sum of finite size and density correction 

5. 

fa 

= 

Aircraft noise frequency 

6. 

aspla 

= 

Actual sonic environment one-third octave spl's on 
aircraft structure 

7. 

sea 

= 

Transducer finite size correction , for aircraft noise 

8. 

mspla 

= 

Measured one-third octave spl for aircraft noise 
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TABLE 8. COMPUTER OUTPUT; S-3A SPECTRUM - LOCATION 1 
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TABLE 9. COMPUTER OUTPUT: S-3A SPECTRUM - LOCATION 2 
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TABLE 10. COMPUTER OUTPUT: S-3A SPECTRUM - LOCATION 3 
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TABLE 11. COMPUTER .OUTPUT: L-1011 SPECTRUM 
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TABLE 12. COMPUTER OUTPUT: V/STOL SPECTRUM 
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TABLE, 13. COMPUTER OUTPUT; SCAT 15-F MODEL TEST - LOCATION 31P-10 


oo to CT> ft 00 «-- 4 o^or*.cM-^oc^oooocM 


m 


d ca 

<rt *o 


-3- LT\ ft tr» O fl 0>CT>00cr>Or^C«4f^C0^cD 

ft ft ft rH iH iH rH ft fH ft *~t iHrHiHfHfHrHpHiHrH 


OOO OOO OOO ftftftftftCstCNrOf^^tO 

tneo ••»••••••••••••••••• 

o T3 OOOOOO OOO ooooooooooo 
(AW I I I I I I I I I I I I I I I I I I I I. 


<0 r^cnLDo^j“Lnr*...ftc7> ftftrvooK\incNi..Jf(Cstoo 


O.C3 
(A “O 
(D w 


^lA-3- ftftCTJ Often (DO^COO^Or^CsJ LAtO 

tAi/\LALrtLA^t/>i/\^ j“-:r^-3'LnLr\tDirt-^-=f.^ 


ru M 


OC CC 

o 

OJ ro 

4J 

M- 
4J O' 

o u 

0) 4; 

4- X 

bO bO 

(A % 

M> (A 

(A (A 


oj a 

JD M- 

*j orjD 

4-» 4^ 4-f C 


■o T? 

f- stfc 

U. (A — 

LW M- 4- (A «<- 





1- 

ft o 

o 


rA o 


o 


4-»*0 

U. 

CM O CM 

00 CD 

OO 

CM CM o 


lA 



< 

-O OO rM O O 

0> LA 


o rA o o J- 


CM 



o: 

• • • • • 

• • 

'• 

« • • « « 


« 



o 

LA LA rA ft O 

LA 

LO 

O ft lA J’ J" 

rA 

O 



a: 

^ cn ft 



o ft cn 

« 




__ 

CM to ^ 



r*«. r>. 



CD 


< 

CM ft 



CM ft 

o> 


■o 



00 


(lO 

t 


E CC 
OTJ 
(A w 


CM 

CM 

rA 

o o> 

CP 

O OO 

CM LA 

LA 


O 

CM 

00 

O 


LP 

o 

r*». 

C7> ft 

J- 

OO CM 

OO 

CD LA 

r«*ft 

rH 


rA 

O 

OO 

LP 

O 

fM 

CM 


ft 

rH 

ft CM 

CM 

rA ^ 

LAh** 

CT> 

ft 


OO 

CM 

00 

LP 

LA 

r*«- 









ft 

ft 

rH 

CM 

CM 

rA 

-0* 

LA 


m-0'..:tK%cstoeo(OLA-:fCMpHh»LAcooo^^-:f^ 

t£>u>^iD(0(OLnLrMA Lnm ^cMCMOrAOOKvrt^ 
iiiiiiiiiiiiiii pHfipH 


rMCMrO-^tAr^OOftCsJmUMOOCNOlAftftOrH 


OOOOOOOftftftrHCM^^tOOKW^OO 

ft iH ft CM CM 





c 

O' 


c 


II 




f 

(/I 


•— 











c 

eo 


o 






o 

LA 


o 


o> o 




•— 

LA 


o 

o 

00 ft 



CO 



OO 

o 

OO O 

^ o o 

ft 

o 

ft 

u 

O 00 o> 

lA 

o 

oi o 

ft iH -If LA 

tA O 

CM 

01 

• • • 

• 

• 

• • 

« • • • 

• 

• 

• 

L. 

o ^ rA 

lA O 

CM 

3 

0 

1 

0 

o> 

^ ft 

o 

L. 

^ CD 




LA 

cn 

(A • 


o 

LA rA 





00 

<7» 


CJ 
















iH 


> 









tJ 








(0 

»— 





fO tu 

> 

4-f 


V) 



X 

<0 

® O E - 

c 

VI 

•o 

c 

U </i u 




to 0) 

o 

— c 


Q) 

*J *J Q. 

2 

E 

u*o 

(0 U bO > 

u 

*o u_ 

4D 

o 


000>lOCO«MrMCNfOOrMO>4t^ft.^COOlOfAlAO CO 

^03 •••••••••••••«•••••• *0 

O.TJ miOLACMCMOfifioc7>oors*r««.r«%cMrAftrM(Ooo 

(Aw (O 

fM rM rM i-M i-M rM rM rM ■ rM f^ rM rM 


e N 

4- X 


OOO 

irt o o 

CM O O 
ft ft CM 


OOO 

O LA o 

LA rH O 
CM ro 


OOO 
OOO 
O K> O 
LA LD CO 


OO OOO 

OO OOO 
OLA O o O 
OCN O O LA 
ft ft ft CM CM 


O O O O O 

O O O o o 

LA O O O O 
ft O O rA O 
rA LA LO 00 


o 

o 

o 

o 

o 


92 


0ASPL= 166 



TABLE Ik. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 31P-12 
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TABLE 15 . COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 32P-10 

MODEL AIRCRAFT 
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TABLE 16 . COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 32P-12 

MODEL AIRCRAFT 
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TABLE IT. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 33P-10 
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TABLE 18. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 33P-12 
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TABLE 19 . COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 20P-10 
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TABLE 20. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 20P-12 
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TABLE 21. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 92P-10 
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TABLE 22. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 92P-12 
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TABLE 23. COMPUTER OUTPUT: SCAT 15 -F MODEL TEST - LOCATION 36 P-IO 

MODEL AIRCRAFT 
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TABLE 24. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 36P-12 
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TABLE 25. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 37P-10 
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TABLE 26. COMPUTER OUTPUT:' SCAT 15-F MODEL TEST - LOCATION 3TP-12 
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TABLE 2T. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION 46P-10 
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TABLE 28. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION U6P-12 
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TABLE 29 . COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION LTP-10 

MODEL AIRCRAFT 
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TABLE 30. COMPUTER OUTPUT: SCAT 15-F MODEL TEST - LOCATION kjP-12 

MODEL AIRCRAFT 
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